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' We will analyze the itinerant model for ferromagnetism with both single-site and two-site electron 

correlations. We will include band degeneration into the model. This will allow us to consider the 
■ on-site exchange interactions in the Hamiltonian. The modified Hartree-Fock approximation for the 

' two-site interactions will be used. This approximation will give the relative one spin band broadening 

with respect to the other, in addition to the shift in position of majority and minority spin bands. 
Next, we will use the coherent potential approximation technique (CPA) with two self-energies. One 
will describe the on-site correlation and the second one the inter-site correlation. We will separate 
from these self-energies, the linear terms arising from Hartree-Fock approximation and the higher 
order terms. The on-site linear term will contribute only to the effective molecular field. The 
inter-site linear term will contribute to the effective molecular field and to the inter-site correlation 
factor. The Green function technique and CPA decoupling will allow for the change in the shape of 
C/2 ' spin bands, which has been described by the correlation factors and which will decrease the kinetic 

energy of the system. The gain in kinetic energy due to the on-site and inter-site correlation factors 
will drive the ferromagnetism and significantly reduce the effective Hartree field necessary to create 
this ordering. 

a 

o . 

1. Introduction 

The basic model for magnetic ordering of itinerant electrons in solids is the Hubbard model [1]. In the mean- field 
approximation, the Hubbard model leads to the well-known Stoner model for magnetism [2,3]. The Coulomb constant 
U coming out of the Stoner condition for creating ferromagnetism is large, i.e. of the order of bandwidth. On the 
one hand it can be justified by the existing strong Coulomb interaction, but on the other hand for such a strong 
interaction, one can not use mean-field approximation. This has prompted attempts to treat the problem within the 
higher order perturbation theory. From the many approaches considered by others, we mention only Hubbard I and 
conventional CPA [4] , which like many other approaches have failed to bring any type of ferromagnetic ordering [5] . 
"j—s ■ They did not produce the spin-dependent band shift necessary for a ferromagnetic ordering which is why new versions 
of the conventional CPA are still being created. One of these new attempts is a self-consistent moment method called 
"spectral density approach" (SDA) [6]. The magnetic phase diagrams calculated by this method are more realistic 
and the values of the Curie temperatures also make sense. The main advantage of SDA is obtaining a spin-dependent 
band shift. The main shortcoming of this method is the lack of quaziparticle damping. Therefore, Nolting and co- 
workers have proposed a combination of SDA and CPA called "modified alloy analogy" (MAA) [7]. Nonetheless, in 
O ■ our knowledge there is no reliable solution of the magnetic ordering problem in the Hubbard model, and still to this 
day the problem has remained unsolved. More recently, the new dynamical mean-field theory (DMFT) [8] has been 
developed for a direct computational simulation of systems with correlated electrons on a crystal lattice. This method 
has an exact solution in the non-trivial limit of an infinite coordination number [9] . The results have been obtained by 
using direct quantum Monte-Carlo (QMC) simulation and the mean-field Green function theory. Use of this DMFT 
method, has introduced a significant progress in the theory of ferromagnetism. The results (see e.g. [10,11]) show the 
existence of ferromagnetism but at much lower temperatures than those coming from the Hartree-Fock approximation. 
Such results would remove the problem known as a "magnetic paradox", i.e. the Curie temperature coming from 
interaction constant U -fitted to obtain the proper magnetic moment at T = OK- as being much too high. In this paper 
we will use a more conventional approach of many body theory. We included in the CPA method both the on-site 
and inter-site correlations. We will show, despite the use of a traditional approach, that qualitatively new results can 
be calculated which can bring the constant of the mean field creating magnetism to almost zero. Extension of the 
Hubbard model is the model which, in addition to the on-site Coulomb repulsion U includes the on-site exchange 
interaction J„ = (iX, iv \ l/r \ (for A ^ z/, where A, v are the sub-band indices) and the inter-site exchange 
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interaction J = {iX,jv 1/r jv,iX) [12]. To obtain a magnetic state, many of the previous papers on metallic 
ferromagnetism have utilized the on-site exchange interaction (see e.g. [12,13]). Hirsch [14-18] and others (see e.g. 
[19-21]) included the nearest-neighbor inter-site interactions without the band degeneracy; | 1/r \ k,l), where 
(k,l) — and i, j are the nearest-neighbor lattice sites. Their solutions of the 3 dimensional Hubbard model did 
not yield the exchange field necessary for ordering which would have been small enough to be justified physically. 
In our model, we introduce into Hamiltonian the on-site Coulomb repulsion U, single-site two sub-band interactions 
(i\,jv\l/r\k'd, Iw) (i = j = k = I and for the sub-band indices we have the restriction (tn, i9) = (A, v) and A ^ v), and 
also the inter-site nearest-neighbor interactions (i, j \ 1/r \ k,l) (for (k,l) = and i ^ j). Analyzing the inter-site 
correlations, we deal with the probabilities defined for the product of four operators (the two operator product is 
replaced by its stochastic value). We have interpreted stochastic value of the two operators product as the product of 
probabilities for electron transfer between two atoms in the presence of another electron. This interpretation makes 
the connection between Hirsch's [14] average bond occupation for spin a; I a , and the standard probabilities used in 
the CPA method [4]. This new link allows us to expand the model to the more realistic cases of the weak intcr- 
site interaction being accompanied by the single-site Coulomb interaction of any strength. The spin band narrowing 
(broadening) is arising from the two-site interactions (see [14,22]) and also from the on-site interaction (see [5,23]). On 
the basis of general understanding, how the correlation affects the density of states (DOS) (see e.g. Fukuyama [23]), 
we can interpret the part of the band broadening which is coming from the two-site interactions, as being expressed 
by the inter-site correlation factor Kij, even in the lowest-first order approximation. 

The paper is organized as follows. In Section 2, we have put forward the model Hamiltonian and developed the 
CPA formalism to treat the on-site and inter-site Coulomb correlation at the same time. In Section 3, we have 
set up the model for ferromagnetism, which includes the on-site Coulomb correlation, inter-site correlation and the 
assisted hopping correlation. Numerical examples are presented in Section 4 based on the semi-elliptic DOS for the 
weak inter-site interactions, in the presence of on-site Coulomb interaction of any strength, with and without hopping 
correlation. On the basis of these results, the conclusions regarding the appearance of magnetic ordering with growing 
occupation of the band are drawn in Section 5. 



2. Hamiltonian and the Coherent Potential Approximation 

The Hamiltonian for one degenerate band can be written in the form given by Hubbard [24] 

H = - Uj (c+ Xa c 3 „ a + h.c) - fi a Y + (i\jv\l/r\kti, M c^ CT c+ CT ,c to(T /c M(T , (1) 

i,\,<7 ijkl, 

a Xuvoi) 

i 

(7,0- 

where Uj is the nearest neighbors hopping integral, fio is the chemical potential, cf Xa . (ci\ a ) creates (destroys) an 
electron of spin a in a Wannier orbital A on the i-th lattice site, the indices A, v, w, 1? numerate the sub-bands in the 
degenerated single band. Taking into account in Hamiltonian (1) only single-site i = j = k = I and two-site interac- 
tions ((k, I) = as well as single sub-band (A = v = w = §) and two sub-band interactions ((w, = (A, v)) we 
have obtained and retained the following matrix elements; 

• single-site, single sub-band interaction 

U = (iX,iX\l/r\iX,i\), (2a) 

• single-site (subscript " in" ) , two sub-band interactions (for A^v) 

Vin = (iX,iv\l/r\iX,iu) , J in = (iX, iv\l/r\iv, iX) , J' in = (iX, iX\l/r\iv, iv) , (2b) 

• two-site interactions 

V = (iX,jv\l/r\iX,jv) , J = {iX, jv\l/r\jv, iX) , 

J' = {iX,iX\l/r\jv,jv) , At - (iX,iX\l/r\jv,iX) . (2c) 
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The mean- field approximation was applied to intra-atomic weak interactions Ji n , J' in , Vj„ and to the assisted hopping 
interaction, Aio- After this operation, this part of the interaction together with the kinetic energy of Hamiltonian (1) 
took on the following form 

H = -J2 % {4 a c ja + h.c.) - Vofiia + J2 M "^a, (3) 

ij ia i<r 

a 

where the molecular (mean) field is given by 

M a = In-" + 2zAtI- (T (m), (4) 
with / being the effective on-site exchange interaction constant 

/ = (d - 1) {J in + J' in + V m ) , (5) 

where d is the number of sub-bands, and I- a (m) the average bond occupation for spin — a and magnetization m 

I- a (m) = (cl_ a Cj- a ) , (6) 

where z is the number of nearest-neighbors. The spin-dependent hopping integral tfj is expressed by 

t^=tij- M{ni- a + nj- a ), At = d-At . (7) 

The approximation given by Eq. (7), changing the bandwidth, was proposed by Micnas [22] and Hirsch [14,15]. 
After the Fourier transform we obtain 

H Q = J2(4-»)rika, (8) 

k<7 

where 

el = e k (l - 2^n" CT ) - a^m, (9) 

s k = -ts k , s k = e ik[Ri ~ Ri) i M = Mo-^n (10) 

<i,j> 

and t = Uj is the nearest-neighbor hopping integral. 
The modified molecular field is 

M M° 

M a = —n - a—^m, 11 
2 2 v ; 



where M ( m ) have the following forms: 



2zAi , N 

M = /+ 1 , (12) 

n 



M^ = I-a—(I- a (m)-I ), (13) 
m 



where 



I Q = lim J±cr(ra). 



We will leave now, for the time being, the single-site Coulomb repulsion and only consider, in the Hamiltonian (1). 
the inter-site terms and small intra-atomic terms; Ji n , J' in , Vi n , treated in Hartree-Fock approximation. The following 
form was arrived at 
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ka <ij>a 

where the new inter-site constants J, J', V for the degenerate band are given by dJ a , g?Jq, dV respectively, where d 
is the number of degenerated orbitals in the band. The physical meaning of these interactions is the following; J is 
the exchange interaction, J' the pair-hopping interaction, V the density-density interaction. 

The CPA idea was applied to the Hamiltonian (14) and the stochastic value replaced each operator product in the 

square bracket. For example, the operators product; Cj G Ci a was replaced by the stochastic value c^ a Ci a . It took the 
values lor 0. The probability of value 1 is the probability of the electron with spin a hopping from the i to the j 
lattice site which is given by the product of probabilities that there is an electron with spin a on the i site and that 
the j site is empty; n? (1 —n?). This probability is called Pf below. 
Introducing the inter-site self-energy Sf 2 it can be written that 

H = J2 E° k n ka + ]T (ei - £?, 2 ) c+cjr, (15) 

a 

where the effective dispersion relation (see Eqs (9) and (10)) is given by 

E% =e%-n + flfcEf, = e k (l-^- 2^n-A - o M ™ + zJ m - ,i. (16) 



According to what was said above and from Eq. (14) the stochastic potential £j is equal to 



It will take on one of the following values 

(e 1 = J-V Pf = I a (m) = n? (1 - n?) 

£i = I e 2 = J+J' ; with probabilities; P 2 °" = I-„{m) = nj a (l - n~ a ) . (17) 
{ e 3 = Pg = l-P?- P% 

Pf and Pf above arc equal to the quantities I a , l_ a in Rcf. [14], if we assume the ferromagnetic order, i.e. n\ = 
rij = n a . Once again, e.g. the probability Pf is the probability of electron with spin — a hopping from the j to the i 
lattice site in the presence of another electron with spin a (hopping from the j to the i site) detected by the operators 
product c+ Cj CT , which is at the end of Eq. (14). 

With these probabilities and potentials one can write the following equation for the self-energy 2 i which describes 
the inter-site interactions only (Hamiltonian (15)) 

fr v^c. -* (18) 

where Pf and Si are given by Eq. (17) and the Slater-Koster function has the following form 

°-w4£ .-^U^ (19> 

In the first approximation, it was obtained from Eq. (18) that 

S? i2 ^^ i2 = P 1 ^i+P 2 CT e 2 . (20) 

The effective dispersion relation in the first order approximation is given on the base of Eq. (16) and (20) by the 
relation 



E k — £ k + s fc e l,2 - A* 

= £k 



D D D ' 



-a — m-fj,, (21) 
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which is identical to Eq. (4) in Ref. [15] with the exception of the At term included here. One can see that this 
result is only the first-order approximation of the self-energy Sf 2 m J/D, J' / D and V/D, where D — zt is the 
half-bandwidth and z is the number of nearest-neighbors. Equation (18) allows for the higher order calculations of 
the inter-site interactions through the inter-site self-energy £f 2 . 

It is worth mentioning here that in the first order approximation (with respect to the differences (ej — e{) / D) the 
self-energy EJ 2 in Eq. (20) depended on the spin index a, but as is known today (see Refs [5,23,25]), the mean-field 
approximation tends to overestimate the effects of ordering. When summed up, the infinite CPA series, coming from 
Eq. (18) and the self-energy dependence on spin index disappeared, i.e. Sf 2 = ^12 = ^1,2- As a consequence, the 
full solution of Eq. (18) has given the ground state which remained paramagnetic until an additional molecular field 
coming from the Hartrcc-Fock approximation (used for some single- and two-site interactions) was introduced into the 
Hamiltonian (1). This molecular field brought the energy dependence on u through the field cr^rpm (see Eq. (13)). 

To include effect of on-site Coulomb repulsion, in addition to the two-site interactions, it was now assumed that 
the effective stochastic atomic energy Ef at the i-th site is given by (see Ref. [23]) 

£?? = E& + e<) (22) 

where £i is given by Eq. (17) and the probabilities of energies Ef are also given by Eq. (17), but now the Eq. (18) is 
replaced by 

Jt, Ef-YX 

YP? r J —s = 0, (23) 



with 



and with the single-site self energy given by 

_ Un~" 

where the zero energy is defined at the atomic level. 

It is worth noting here that because of the additional factor Sk in Eq. (24), there was a modification to the standard 
method of calculating the DOS. When we used the relation Sk = —£k/t, and the relation (9) for we were be able to 
write 



1 

N N " 



te-e k (l- 2fn-° - Ha) + a^+llm - E£ + „ 



1 / w , MZ.+zJ , \ 



where the unperturbed Slater-Koster function is given by F (e) = j?Yl j"rj- an d 

k 

^ = l_M_ 2 ^ n - (27) 

In our derivation we made the assumption that both self-energies were ^-independent, while they can be energy- 
dependent. 

3. Ferromagnetism (F) 

At nonzero temperature the number of electrons with spin a is given by 

00 

n a = - j f (e) ^ImF CT (e) de, (28) 
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where /(e) is the Fermi function 



f (e) = — a \ ( I w \\ (29) 

and 

'■w-sErns^w <30 » 

The energy e in above expressions is from now on equal to e + /i + <rm{M^ n + zJ)/2. 

To obtain the criterion for ferromagnetic state we differentiate Eq. (28) with respect to to. As a result, we get 
two terms. The first one is the change of the DOS with rising to. It depends on the self-energies Ef 2 and Egr and 
contributes to the correlation factor K, see Eq. (32) below. The second term is the spin-dependent band shift, which 
contributes merely to the Hartree field. In effect we obtain the following formula 

oo 

l = K+[I cr + zJ + zAt(l-n)} J p° m=a (e)P T (e)de (31) 

— OO 

where 

CO 

p a (e) = —TmF°(e), P T (e) = f (e) e^~^(3, K = --Im [ f( £ )^-^-de (32) 

— oo 

and J cr is given by Eq. (5). In the zero temperature limit the function Pt (e) becomes the Dirac's delta function and 
we have 

l = K + p{e F )[I cr + zJ + zAt(l-n)}, (33) 

where 

p(e F ) = -ilm F^)| m=0 , = -|lm J ^j^ds (34) 

— £> 

and Sf is the Fermi energy calculated from Eq. (28) in the zero temperature limit. 
Using relation 

dF a dF 17 dZf T dF a d^ 2 dF° da° , N 
= \L _i j (35) 

dm <9X^ dm <9E? 2 dm da" dm 1 v 1 

where a a — 2^rT a , we obtain the correlation factor as the sum of the on-site, inter-site and assisted hopping 
correlation factors 

K = K u + K ij +K At , (36) 

where 



2 e jdF°d?% 2 J dF° d^ ia 

Ku = -n lm J d^^^ = -^ Im 7 gE ?2 9to & ' 

-D —D 



2 T T 9F a <9a CT , 

-D 

To establish the link with the existing literature, it can be written on the basis of Ref. [23] for the single-site correlation 
factor in the case of strong correlation (U 3> D) that 
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Kjj = Imln 



For the two-site correlation factor, using Eq. (30) in Eq. (37), we obtain 



Kij = Im 



i asf,2 (f) 

6^ 9m 



cfe, 



(38) 



(39) 



where 

Later on we used the explicit form of F a (e, m), obtained after the elimination of self-energy Although this 
elimination is only possible in some specific cases, e.g. the semi-elliptic DOS, there is a great advantage of the 
simplicity in this case since there is no need to develop a formula for the correlation factors any further than has been 
formulated in Eq. (34). 

For the assisted hopping correlation factor, using Eq. (30) in Eq. (37), we obtain 



K At 



2 At 



Im 



1 



de, 



(40) 



-D 



where Ffi was defined earlier after Eq. (39). 

In the lowest order Hartree-Fock approximation; Sf 2 
of Eqs (39) and (40) we obtain 



e l,2> n U 



Un- a , Eg- = and from Eq. (33) with the help 



1 = p{e F ) 



I cr + zJ+U + — 



2e F 
T 



ut l,2 

dm 



-At 



(41) 



The factor I cr + zJ + U above comes from the Hartree-Fock approximation, the extra term in the curl bracket; 
(2e F /t) ((0ef )2 /0m) - At), comes from the inter-site and assisted hopping correlation factors. These factors collect 
contributions to the magnetic criterion coming from the change of DOS over the whole energy interval, as opposed to 
the contribution coming only at the Fermi energy, as in the case of the classic Hartree-Fock approximation (the first 
Stoner type term in Eq. (41)). 

It is worthwhile to note here that in the first order approximation; T,fj w U t nT a and the on-site correlation factor 



Kjj is equal to zero. Without this on-site correlation factor (Kjj — 0), together with SJ 2 
we obtained from Eq. (41), for the constant DOS, the Hirsch's result (sec Rcf. [15]) 

2 [1 - (I?r + U)/D] 



-1,2) 



l + a-(3a-l)(l-n) 



2 ' 



as well as At = 0, 



(42) 



where aj = = § and j = = §. 

From Eq. (41), without the correlation factors, the Stoner criterion of magnetism was reached in the following 
well-known form 



1 = p° (e F ) (I cr + U) 



(43) 



where p a (e F ) is the unperturbed DOS on the Fermi level. 

Numerical results in the next section were calculated on the basis of equation (33), which was obtained in the zero 
temperature limit of Eq. (31). The correlation factors in this limit are defined by Eqs (37). 

4. Numerical examples 



This new model will be illustrated by showing the dependence of the critical on-site field versus the carrier concen- 
tration; I cr (n). It must be remembered that the relation between this field and the total molecular field M computed 
here is as follows I = M — zAt(l — n) — zJ. Only the first order approximation was considered for the two-site 
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self-energy 2 = e\ 2 . The single-site self-energy was solved using the standard CPA equation (see Eq. (25)). 
The perturbed Slater-Koster function has the form 



where F (e) is the unperturbed Slater-Koster function, for which we assume the semi-elliptic form 



Fo(s) = - 



£ 

D 



DJ 



(44) 



(45) 



"1.2 
/ 



2fn~ 



and b a = 1 

To illustrate the new results we used Eq. (33) to calculate the dependence of critical (minimal) value of I cr on the 
carrier concentration n. The inter-site correlation factors Ky and K^t are calculated from Eq. (39) and Eq. (40) 
with the help of function F a defined by Eq. (44). Inter-site interactions e\ and £2, parameter p and the asymmetry 
parameter a are defined as follows 



ei = ( J — V) = ap — = apt and e% = J + J' = p — = pt. 

z z 



(46) 



We assume that the inter-site exchange and pair hopping interactions are equal J — J' (see Ref. [19]), then J = 
pD/2z =pt/2. 

In all the figures presented below, the dependence I cr (n) is expressed in units of half-bandwidth D. 




FIG. 1 



FIG. 1. Dependence of critical field on carrier concentration I cr (n), for the single-site correlation alone 
(U 0, J = J' = V = At = 0), for different strengths of U; U = 5D is the solid line, U = D the dashed line, U = 0.5D the 
dotted line. For comparison we show the I cr (n) dependence without on-site Coulomb correlation - dot-dashed line and for the 
strong single-site correlation; U/D = 00 is the double dot-dashed line. 



The function I cr (n) for the single-site correlation (U 7^ 0) alone ( J = J' = V = At = 0) is shown in Fig. 1. The 
solid line corresponds to U — 5D, the dashed line is for U — D, and the dotted line is for U = 0.51?. For comparison, 
it is presented in the same figure the dependence of I cr (n) without on-site Coulomb correlation (the dot-dashed line 
for Kjj — 0, which is the standard Stoner criterion for magnetism) and for the infinitely strong single-site correlation; 
U/D = 00 (the double dot-dashed line). 

Analyzing the I cr (n) curves presented in Fig. 1 one can see that the on-site correlation Kjj decreases the minimal 
value of atomic field I cr necessary to obtain the ferromagnetic state. The dominant reduction takes place at the 
half-filling, already at relatively small U « 3D. Further increase of U does not influence the ferromagnetic state 
much. 

In Figs 2-4, we have presented the dependence of critical on-site interaction I cr (n) in the presence of inter-site 
correlations. Figs 2 and 3 do not have hopping interaction (At = 0). With the symmetric initial DOS (semi-elliptic), 
the curves for these interactions remain symmetric with respect to n = 1. This is why the results in Figs 2 and 3 are 
drawn only for < n < 1. 
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FIG. 2. Dependence of critical field on carrier concentration I cr {n), at U = 3-D, a = —1.5 and At = for different strengths 
of inter-site correlation (different parameter p); p — 1 is the solid line, p = 0.5 the dashed line, p = 0.25 the dotted line. For 
comparison we show the 7 cr (n) dependence for p = the double dot-dashed line, and for both inter-site correlation p = and 
on-site correlation Kjj — is the dot-dashed line. 

Fig. 2 shows the influence of parameter p = z(J + J')/D = (J + J')/t on the value of on-site critical field I cr . 
The solid line is for p = 1, the dashed line for p — 0.5, and the dotted line for p = 0.25. The Coulomb interaction; 
U = 3D and the asymmetry parameter of inter-site interactions; a = = -j-fjr — —1.5. We have shown in the same 
figure for comparison the dependence of I cr (n) for the standard Stoner criterion of ferromagnetism (the dot-dashed 
line) and I cr (n) for single-site correlation only; U = 3D and p = (double-dot-dashed line). Analyzing the presented 
curves, it is possible to see that the inter-site correlations significantly decrease the value of the intra-site critical field 
necessary for magnetism I cr . Equivalently the decrease of parameter p increases the minimal value of atomic field 
I cr . The decrease of p means the decrease of inter-site interactions which causes the decrease of inter-site correlation 
factor Kij and the increase of DOS on the Fermi level p(sf) especially for n close to zero or two. 

3 
2.5 
2 

'U 1.5 
1 

0.5 

°0 0.2 0.4 0.6 0.8 1 

n 

FIG. 3 

FIG. 3. Dependence of critical field on carrier concentration I cr (n), at U = 3D, p = 1 and At = for different asymmetry 
parameters a; a = 0.5 is the dotted line; a — —0.5 the dashed line; a — —1.5 the solid line. For comparison we show the I cr (n) 
dependence for p = is the double dot-dashed line, and for both inter-site correlation p — and on-site correlation Ku — is 
the dot-dashed line. 



The role of factor a in minimising the critical field for magnetism is shown in Fig 3. The solid line is for a = —1.5, 
dashed line for a = —0.5 and the dotted line for a = 0.5. All the curves have p = 1 and U = 3D. We have shown 
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in the same figure for comparison the dependence of I cr (n) for the standard Stoner criterion of ferromagnetism (the 
dot-dashed line), and I cr (n) for the single site-correlation only; U = 3D and p — (the double-dot-dashed line). 
Analyzing I cr (n) curves for different parameters a one can see that the decrease of a causes the increase of I cr for 
concentrations close to half-filling but it decreases I cr at small concentrations and concentrations close to full-filling. 
This effect is the result of change in two factors; inter-site correlation factor Kij and the DOS on the Fermi level p(ef)- 
The value of p (ep) increases with a, but the value of inter-site correlation factor decreases with increasing a. The 
sign of parameter a depends on the difference between interactions J and V. For typical 3-d metal we have J < V 
[1,21] and parameter a is negative. However, when J > V [15] or V < (which can take place in superconducting 
cuprates) the sign of a will be positive. 

As can be seen from Figs 2 and 3, at some values of parameters p and a, we have obtained the ferromagnetic state 
already at zero value of intra-site field. For the optimal values of p = 1 and a — —1.5 the ferromagnetic state exists 
at concentrations n — 0.09 — 1.91 and also at n = 1.78 — 1.91, the latter is more interesting from the point of view of 
ferromagnetic elements. 
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FIG. 4. Dependence of the critical field on carrier concentration I cr (n), for a — —1.5 and p = 1 in the case of correlation 
U = 3D for different values of parameter At; At = 0.2t is the solid line, At = O.li the dashed line, At = the dotted line. 
For comparison we show the I cr (n) dependence for p — 0, At = 0.2t is the double dot-dashed line, and for both inter-site 
correlation p = and on-site correlation Ku — is the dot-dashed line. 



In Fig. 4 we show the critical value of the internal exchange field I cr in function of carrier concentration at a = —1.5, 
p = 1 and for different values of the hopping interaction At in the case of single-site Coulomb correlation U = 3D; 
At = is the solid line, At = O.lt the dashed line, and At = 0.2t the dotted line. We show in the same figure for 
comparison the dependence I cr {n) for p = 0,At = 0.2t is the double dot-dashed line and without both inter-site and 
on-site correlation [K — 0) the dot-dashed curve. The increase of the hopping interaction At causes for all electron 
occupations n to decrease in minimal value of the internal exchange field creating ferromagnetism I cr . However, a 
closer look into Fig. 4 reveals that the largest decrease takes place for carrier concentrations below half-filling. As 
a result, the hopping interaction helps ferromagnetism more in the beginning of the 3-d row of elements than at the 
end. The increase of the ratio; I/D along the 3-d row, as well as the increase of other interaction constants along this 
row (see Ref. [26]), can offset this effect. At half filling, we have discontinuity of I cr {n). Coming from the left and 
right side with n one has 

p(e F ) \ t J 1 h 

with the upper and lower sign respectively. The discontinuity exists only in the split band limit, which takes place 
for our U (our U = 3D). It is absent in the paper [14], since the Coulomb correlation was treated there in the 
Hartree-Fock approximation which can only shift the spin bands and is unable to change their shapes and to split 
them into sub-bands when U is high enough. 
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5. Conclusions 



As can be seen from Fig. 1 the single-site electron correlation, Kjj, seriously helps ferromagnctism, especially for 
concentrations of n close to half-filling. For the 3-d group of elements, the enhancement will not be the same for small 
and large n as we have to take into account the fact that the ratio of I/D existing in different elements growths with 
n. This increase will increase the chances of ferromagnctism towards the end of the 3d-group of elements. 

As was already stated, the strong single-site Coulomb correlation U > D favours the ferromagnetism (but rather 
antiferromagnetism) at half-filled band (see Fig. 1). Comparing this result with the 3-d group of elements, it is possible 
that this situation would correspond to antiferromagnetic elements like Cr, and Mn, which have concentrations of 
electrons close to half-filling. The cohesion energy,£? co h, strongly decreases for these elements which gives evidence 
of the decreasing bandwidth 2D to a small value. Such a decrease of the bandwidth follows from our Eq. (27) with 
increasing interaction constants. In the first order approximation it can be written as; D E g = D[\ — p(l + a)I ], where 
Io has maximum at half-filling. The small value of effective bandwidth would justify using models with U » D for 
these two elements. 

In extending description of ferromagnetism by the CPA method we included; 

(i) spin dependent band shift coming from the intra- and inter-site exchange interactions (I and J, V) and assisted 
hopping interaction At, 

(ii) spin dependent change of the band width depending on electron concentrations and all the interactions. Nar- 
rowing of the band increases DOS on the Fermi level p(sf), 

(iii) band shape changes through the on-site (Kir) and inter-site correlation factors [Kij and K\ t )- 

Figs 2-4 show the role of inter-site correlations for magnetism. Comparing the curves with different parameter 
p = z(J + J')/D in Fig. 2 one can see that the inter-site correlation additionally decreases the intra-atomic field 
necessary to create the ferromagnetism. The value of on-site Coulomb repulsion was fixed at the realistic level of 
U = 3D. The effect is particularly strong for large but still realistic p — 1. One can see that the decrease of intra- 
atomic field (symmetrical in n = 1) is even larger than the one coming from the single-site correlation Kjj and it 
takes the place for values of n corresponding better to ferromagnetic elements. The influence of different inter-site 
interactions on ferromagnetism is given by the parameter a = = j-fjr- 

It is strongest at both ends of the band and for a as small as possible, we assume a = —1.5, see Fig 3. Such a value 
is possible for transition metals, since for these elements values of inter-site charge-charge interaction V are positive 
and larger than the inter-site exchange interaction J. For a = 1, i.e. J + J' = J — V the inter-site interactions do not 
change the effective bandwidth (see Eq. (21)), hence the inter-site correlation factor Kij = and the only influence 
of inter-site interactions on ferromagnctism is by increasing the total molecular field M. 

To sum up, we can say that the inter-site correlation in co-operation with single-site correlation enhances ferromag- 
netism preferably at the end of the band, or speaking about 3-d row, at the end of this row. 

Fig. 4 shows that the hopping interaction At also enhances magnetism, but this effect is stronger for smaller 
concentrations of n than for larger concentrations. This result seems to be in contradiction to the experimental 
evidence on the 3-d row of elements, but the increase of the ratio I/D and of other interactions along the 3-d row of 
elements can offset this effect. 

In general the CPA method used for the Hubbard model, without any additional field of on-site or inter-site origin, 
does not yield the ferromagnetic state. It comes out of the fact that the spin densities are not shifted or twisted with 
respect to each other (e) = pi (e). For comparison, the Spectral Density Approach (SDA) [6], already mentioned 
above, gives the ferromagnetic state for Hubbard model at U > D. The main advantage of the SDA is obtaining 
spin band shifts due to 'higher hopping t%j correlations'. The drawback of this method is that the self-energy is real 
which neglects quasiparticle damping, that takes place in these materials [7]. Comparing the SDA method with our 
model, one notices that the so called "higher" correlation function (defined as B- a in [7]) is the interaction very 
similar to the assisted hopping interaction At used in our paper. Using moments method, introduced by Harris and 
Lange [27] and more recently by Herrmann and Nolting [7], the CPA technique was applied to two modified atomic 
levels with modified probabilities, by the spin depended shifted atomic level -B_ CT . This approach was called Modified 
Alloy Analogy (MAA). As a result, Herrmann and Nolting obtained within MAA, the ferromagnetic state for some 
concentrations with the maximum of magnetic moment around n « 0.7. Using the SDA method, the same authors 
obtained ferromagnetism for concentrations n > 0.55. Using our model, with only the assisted hopping interaction 
At, we obtained very small values of the on-site exchange interaction / necessary for ferromagnctism at n « 0.7 (see 
double dot-dashed line in Fig. 4). This result is very similar to the MAA result mentioned above. Introducing into 
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consideration the other inter-site interactions, we got ferromagnctism for smaller concentrations and even at zero 
values of the on-site exchange interaction I (see Fig. 4). 

The ferromagnetic state in the standard Hubbard model can be obtained for the infinite coordination number 
(Z — > oo) using dynamical mean-field theory (DMFT). This method, if used for symmetrical DOS (like the one 
used by us), needs very large values of on-site Coulomb correlation U. For intermediate U in DMFT method to 
obtain the ferromagnetic state, one has to use strongly asymmetric DOS (see e.g. Ref. [10]). The use of such DOS 
in CPA would also allow for a large reduction of the inter-atomic and intra-atomic field necessary to create the 
ferromagnctism. Vollhardt and co-workers [10,28] pointed out that use of the inter-site exchange and other nearest- 
neighbour interactions, which was done in our paper using a different method, would reduce the value of U necessary 
for fcrromagnetism in DMFT method. 

In general the Hartree-Fock approximation overestimates the ordering. As a result, the same interactions when 
treated in higher order approximation, required much larger interaction constants to obtain the alignment. On the 
other hand, use of the Hartree-Fock approximation in this paper may be justified since, for many electron occupations, 
the magnitude of the field used, is very small as compared to the kinetic energy characterized by the bandwidth. 
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